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Preface 

This  report  studies  the  linearized  collision  operator 
in  the  Boltzmann  equation  with  repulsive  intermolecular 
potentials   V(r)  =  ar 

It  is  shown  that  for  a  >  2   the  collision  operator 

has  a  purely  discrete  spectrum  and  its  eigenfunctions  are 

2 
infinitely  differentlable  Jif.  -functions  which  are  complete 

2 
in  (^  .   The  proof  relies  on  the  formalism  of  pseudo- 
differential  operators. 

The  report  contains  two  parts :    Part  I  deals  with 
the  special  case  of   a  =  4,   the  Maxwell's  molecules, 
while   Part  II  considers  the  general  case  of  inverse-power 
potentials  with   a  >  2  . 
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Part  I 


Maxwell's  Molecules  --  The  Inverse-Fourth  Power  Potential 


1.   Introduction 

The  present  study  is  concerned  with  the  linearized 
Boltzmann  equation  [  5] 

g|=C[f]  (1) 

CCf]    =    I   d^c^expC-c^)    I   dede|g|l(|g|,6)[f|+f -f-^-f]  (2) 

in  the  kinetic  theory  of  gases  with  repulsive  intermolecular 
potentials 

V(r)  =  a^^/r"  .  (3) 

It  is  convenient  for  our  purpose  to  introduce  a  new 
function 

<p{c)    =  f(c)exp(-|c^)  (4) 

and  consider  the  operator 

L<p   =   -exp(-i<;^)C[(t)(c)exp(ic  )]  .  (5) 

"^  2 

which  is  known  to  be  symmetric  and  non-negative. 

We  shall  show  that  for  a  >  2,  L~  is  a  compact  operator 
from  ^  to  X  J  excluding  the  zero-eigenvalue  subspace.  Con- 
sequently the  spectrum  of  L  is  purely  discrete  and  its  eigen- 

2  2 

functions  are  X  forming  a  complete  set  in  i^  . 

In  presenting  these  results,  we  divide  the  study  into 
two  parts.   The  first  part,  which  constitutes  the  present 
paper,  deals  with  the  special  case  of  a  =  4,  commonly  known 


as  Maxwell's  molecular  model  or  Maxwell's   gas,  while 
the  second  part,  given  in  a  companion  paper  [8]  (hereafter 
referred  to  as  Part  II),  considers  the  general  case  of 
inverse-power  potentials  with  a  >  2. 

The  result  for  the  Maxwell's  gas  (a=^)  is  well-known 
from  the  work  of  Wang  Chang  and  Uhlenbeck  P-1] ,  in  which  the 
eigenfunctions  are  explicitly  constructed.   However,  our 
method,  which  is  capable  of  obtaining  the  most  crucial 
theoretical  results  regarding  the  discreteness  of  the  spectrum 
and  the  differentiability  and  decay  of  the  eigenfunctions, 
does  not  depend  on  the  explicit  construction  of  the  eigen- 
functions and  can  therefore  be  carried  over  to  the  general 
case  of  inverse-power  potentials  with  a  >  2. 

A  brief  description  of  our  method  can  be  given  as 
follows.   We  first  formulate  the  operator 

Ac})  =  (N+L)(})  (6) 

in  the  form  of  a  pseudo-differential  operator,  where  N  is  a 
large  positive  number.   Since  this  operator  is  difficult  to 
handle,  we  further  construct  a  new  operator 

Aj^<\>    =  (N+L^)<i)  (7) 

from  (6)  which  has  the  property  that  A(^  =  ^A   if  <i>   has  the 
form 


(})(c)  =  Y^(c,z)iJ;(|c|)  (8) 

where  Y  (c,z)  is  the  £-th  order  spherical  harmonics  of  the 


vector  c  with  z  as  polar  axis.   It  should  be  noted  that  since 
L  Is  Isotropic,  all  elgenf unctions  of  A  must  be  of  the  form 
(8). 

Our  result  then  follows  from  the  proof  that  A.  Is  invertlble 
and  A~   Is  a  compact  operator  from  sC  to  <5C  ,  <3C  being  the 
square-integrable  functions,  and   ||A~  ||  -»■  0   as  Jl  ->■  «>. 

The  basic  features  of  our  method  can  best  be  illustrated 
with  the  special  case  of  Maxwell  molecules,  since  this  special 
case  is  not  overburdened  with  complicated  computations  and 
lengthy  estimates  that  are  inevitable  to  the  general  case  of 
inverse-power  potentials.   At  the  same  time,  the  result  we 
shall  obtain  in  this  paper  for  the  Maxwell  molecules  provides 
the  key  to  the  case  of  general  inverse-power  potentials. 
This  is  so  because  for  large  c  the  collision  operator  with 
an  inverse-power  potential  for  a  7^  4  can  be  approximated  in 
a  certain  sense  by  an  operator  which  is  very  similar  to  the 
collision  operator  for  the  Maxwell's  gas  (i.e.  with  a 
potential  for  a  =  4).   It  is  for  these  reasons  that  we  shall 
first  present  herewith  the  result  for  the  Maxwell's  gas  in 
an  independent  part,  separate  from  the  general  case  In  Part 
II  for  a  >  2. 

In  the  next  section  a  general  discussion  is  given  for 
the  background  of  the  problem  concerning  inverse-power  law 
Boltzmann  equation.   Section  3  provides  some  material  from 
the  theory  of  pseudo-differential  operators  which  is  needed 
for  our  analyses  in  this  Part  I  as  well  as  in  Part  II.   In 
section  4,-  the  collision  operator  for  the  Maxwell's  gas  is 


formulated  as  a  pseudo-differential  operator.   The  operator 
Ap,  called  the  associated  pseudo-differential  operator-  of  A, 
Is  constructed  in  Sec.  5.   The  main  result  of  this  paper  is 
contained  in  Sec.  6  where  we  prove  that  A.  is  Invertlble  and 
Aj,   is  a  compact  operator  from  ,)L  to  X^ ,  and  ||  A.  ||  ^  0  as  £  -»■  <». 

This  shows  that  A    is  compact  from  _)C  toJJtT^,  or    equivalently , 

—1  ?      -P         -2 

L   is  compact  from  Xj_  to  X,  where  X,  ^-^^  the  square-integrable 

functions  orthogonal  to  the  five  eigenfunctions  of  L  with  the 
zero-eigenvalue.   From  this  we  then  conclude  that  the  spectrum 
of  L  is  discrete  and  the  eigenfunction  exist  and  form  a  complete 
orthonormal  set  of  X-   -functions.   Further  computations  show  that 
the  eigenfunctions  are  actually  in  J .      But  this  will  not  be  done 
here,  since  we  shall  present  a  proof  for  this  fact  in  the  com- 
panion paper  [8]  (i.e.  Part  II)  for  the  general  case  of  inverse- 
power  potentials  with  a  >  2. 


2.   The  Background 

For  the  molecular  model  of  hard  spheres,  the  differential 
cross-section  I(|g|,G)  is  finite  and  the  collision  operator 
L())  in  (5)  is  well-understood,  at  least  qualitatively  [  3]  ,  [  6]  . 
For  all  other  practical  physical  models  [ 5] >  such  as  the 
inverse-power  law  model,  the  Lenard-Jones '  model,  the 
Sutherland's  model,  etc.,  the  interparticle  potential  has  an 
infinite  range  and  consequently  the  differential  cross-section 
has  a  non-integrable  singularity  in  9  at  6  =  ' t '^   •   Because  of 
this  singularity,  the  collision  operator  LcJ)  is  a  non-standard 
unbounded  operator  and  very  little  is  known  about  the  nature 
of  this  operator.   The  only  exception  in  this  connection  is 
the  work  of  Wang  Chang  and  Uhlenbeck  [IIJ  in  which  the  eigen- 
functions  of  L(j)  are  explicitly  constructed  for  the  inverse 
fourth-power  potential  model  (the  Maxwell  gas).   But  their 
method  of  explicit  construction  has  not  been  possible  to  be 
generalized  to  other  infinite-ranged  potentials  or  other 
inverse-power  law  potentials  and  a  general  understanding  of 
the  collision  operator  for  infinite-range  potentials  is 
completely  lacking. 

Several  ways  have  been  used  or  suggested  to  avoid  the 
difficulty  caused  by  the  singularity  of  l(|g|,e).   One  is  to 
use  the  angle  cut-off  approximation  [6]  in  which  l(|g|,0) 
is  set  to  be  zero  in  a  small  interval  p-'"'"  9  >  e  near  p-ir  . 
Another  approach  uses  the  radial  cut-off  approximation  [4] 
in  which  the  interparticle  potential  V(r)  is  set  to  vanish 
for  r  >  R  for  a  large  R.   Still  another  approach  [13J  suggests 


the  use  of  quantum  mechanical  (instead  of  classical) 
cross-section  for  I(|gl,9)  which  is  finite.   This  is  a  method 
to  avoid  the  singularity  of  l(|g|,G)  by  quantum  mechanical 
cross-section  for  the  classical  gas  (at  reasonably  high 
temperatures).   It  should  be  distinguished  from  the  quantum 
mechanical  calculation  of  a  low-temperature  gas  for  which 
quantum  mechanical  cross-section  must  be  used. 

In  all  the  above  approaches,  the  differential  cross-section 
is  finite  and  equation  (2)  can  be  separated  into  two  terms  as 

C[f]  =  -v(c)f(c)  +  K[f]  ,  (9) 

K[f]  =  I  K(c,n)f(n)d^n  .  (lo) 

for  which  a  somewhat  standard  theory  can  be  developed. 

However,  in  the  above  approaches,  both  v(c)  and  K(c,ri) 
become  infinite  as  the  parameter  e,  and  1/R  and  h  (the 
Planck  constant)  tend  to  zero.   This  shows  that  we  have  to 
carry  the  two  terms  in  (9)  which  both  tend  to  infinity  while 
on  physical  ground  the  final  results  of  the  theory  must 
remain  independent  of  the  parameters  e,  1/R  and  h  as  they 
become  vanishingly  small.   This  clearly  indicates  that  these 
approaches  cannot  lead  to  a  satisfactory  theory  of  classical 
gases.   Even  from  the  viewpoint  of  numerical  analysis,  the 
form  given  by  (9)  as  the  difference  of  two  extremely  large 
terms  would  not  be  suitable. 

We  then  conclude  from  the  above  consideration  that  to 
develop  a  satisfactory  theory  of  classical  gases,  the 


difficulty  arising  from  the  singularity  of  I(g,e)  cannot  and 
should  not  be  avoided.   The  purpose  of  the  present  study  Is 
to  develop  a  theory  which  Is  capable  of  handling  the  Boltzmann 
collision  operator  L<p   with  I(|g|,G)  having  a  non-lntegrable 
singularity  In  9. 

In  studying  the  collision  operator,  we  are  confronted 
with  the  alternatives  of  working  with   (1)   the  eigenvalue 
problem  for  L(j),   (2)   initial  value  problems  (including  the 
sound  propagation  problem)  of  the  Boltzmann  equation,  or 
(3)   simple  boundary  value  problems  of  the  Boltzmann  equation. 
We  have  decided  to  choose  the  eigenvalue  problem  since  it 
is  the  simplest  mathematically  well-defined  problem  while 
the  remaining  two  types  of  problems  are  more  difficult  and 
therefore  probably  need  some  approximations  which  we  do  not 
desire  to  go  into  at  this  stage.   At  any  rate,  the  eigenvalue 
problem  seems  to  be  worthwhile  to  be  looked  into  and  it  is 
hoped  the  understanding  of  the  eigenvalue  problem  will  be 
useful  to  more  difficult  problems  of  the  Boltzmann  equation. 


3.   Some  Preliminaries 

This  section  gives  a  simple  discussion  of  pseudo- 
differential  operators  (PsDO)  which  will  be  sufficient  for 
our  need  in  this  Part  I  as  well  as  Part  11. 

We  shall  throughout  this  paper  use  c  and  k  to  denote 

three-dimensional  vectors  while  their  scalars  are  denoted  as 

|c|  and  |k|.   We  shall  also  use  the  shorthand  notation 

m,   m„   m^ 
8^  =  8  1  3  2  3  3 


|m|  =  m-,  +  mp  +  m- 


where  c.,  i  =  l,2,3  are  the  three  components  of  c.   Similar 
notation  will  be  used  for  the  k-derivatives . 

An  operator  P  defined  by 


P(})  =  f  e^^*^p(c,k)$(k)dk 


(11) 

is  called  a  PsDO  if  p(c,k),  called  the  symbol  of  the  operator 
P,  satisfies 


3"^3"p(c,k)|   <   a   (l+c2)^(l+k2)^ 


c  k 


mn 


(12) 


for  |m|  >_  0  and  |n|  >_  0  where  a   and  y  may  be  positive  and  (}) 
is  the  Fourier  transform 


({"(k)  =  (2tt)' 


■3 


-1  k  •  c  ,  /  X  , 
e     ())(c)dc  . 


(13) 


In  the  case  of  Boltzmann  collision  operator,  the  indices 
o  and  M  in  (12)  are  positive  and  this  makes  this  class  of 
PsDO  difficult  to  treat  since  differentiation  with  c  and  k 

10 


does  not  bring  down  the  growth  in  c  and  k;  this  is  essentially 
the  difficulty  that  leads  us  to  consider  the  operator  A^  in 
(7)  in  lieu  of  the  operator  A  in  (6).   However,  for  positive 
o  and  y,  the  following  result  can  easily  be  obtained  which  is 
all  we  need  in  dealing  with  the  operator  A  in  (6). 

THEOREM  1.     The  operator  P  in  (11)  with  condition  (12)  maps 
//into  j?f. 

The  proof  follows  from  the  form 


i   e^^'^(l-9^)J{p^(c,k)V^(k)}dk  (13) 


p^(c,k)  =  p(c,k)(l+c2)~J(l+k2)  ^  (14) 

4;(k)  =  (l+k2)^0(k)  (15) 


where  p,(c5k)  and  its  derivatives  are  square-integrable  over 
R-^xR-'.   This  completes  the  proof. 

We  can  show  that  P  mapsjp/  into  ^  if  (12)  only  holds  for 
m=0  and  0  <_  n.  <_  n*,  n*  is  the  nearest  integer  for  n*  >_  o+2. 
On  the  other  hand,  the  operator  A .  in  (7)  leads  to  a  much 
nicer  class  of  PsDO.   In  this  connection,  our  analyses  rely 
on  two  theorems  due  to  Calderon  and  Vaillancourt  [1,2].   One 
of  their  theorems  may  be  stated  as  follows. 


THEOREM  2.     If  p(G,k)  satisfies  (12)  for  0  <  n.  <  3, 

0  5_  m.  ^  3 ,  i  =  l  J  2 , 3  and  with  a  =  y  =  0 ,  then  the  operator  P  in 

(11)  f rom  ca  to^  can  be  extended  to  a  bounded  operator  from 


11 


^2  x,o  ^ .      Furthermore  '^'^^  ^ -novx^   II  P  II  of  P  satisfl 


es 


IIPII   <   b  Max  I  a   I  (l6) 

'  mn ' 

where  the  maximum  is  taken  over  0  <_  m.  ,n.  <_  3  and  the 
constant  b  is  independent  of  p(c,k). 

For  the  proof  of  this  theorem,  we  refer  to  Ref .  1  . 

An  immediate  consequence  of  Theorem  2  is  that  P  is 
compact  for  a  <  0,  y  <  0.   This  may  be  stated  as  follows. 

Theorem  3.     If  p(c,k)  satisfies  (12)  with  0  <_  m.  ,n^  <  3, 

i=l,2,3,  for  o  <  0,  M  <  Oj  then  P  in  (11)  is  a  compact 

2 
operator  on^fi  . 

To  prove  this,  it  suffices  to  consider 


p^(c,k)  =  ij^^(c,k)p(c,k) 


F  (J)  = 
n 


e^^*''p^(c,k)$(k)d3k 


where  ^    (c,k)  is  a  C  -function  having  the  property 

ip  (c,k)  =  1  for  c^+k^  <  n  and  ^    (c,k)  =  0  for  c^+k^  >  2n, 

n=l,2,...  .   Then  P  are  compact  operators  and 


II  P  -  P  II  ^  0 
n 

as  n  ->•  °o,  according  to  Theorem  2.   Therefore  P  is  compact. 

The  other  theorem  of  Calderon  and  Vaillancourt  that 
we  have  to  use  concerns  PsDO  with  multiple  symbols .   For 
our  purpose,  it  is  sufficient  to  Introduce  these  operators 
as 


12 


Pc|)  =  I  I  e^^'^''"'^^p(c,n,k)4.(n)d3nd3k       (17) 


where  p(cjr),k),  called  a  multiple  symbol,  satisfies 

|9j9j|9|Jp(c,n,k)|   <   p^d+k^)-^^  ,    y  >  0     (18) 

for  all  n.,m.,£.  such  that  0<n.  <4,   0<|ml<8,  and 
0  ^  |Ji|  <  8.   It  is  not  difficult  to  see  that  P,  defined  by 
(17)  and  (18),  maps  ;zf  into ac^  .   In  fact,  we  can  get  an  extra 

(l+n  )    factor  in  the  integrand  by  using  (l+n  )  ()>  e  _ji} ,  an 

2  -2 
extra  factor  (1+k  )    at  the  expense  of  differentiation  with 

2  —1 
respect  to  n  and  an  extra  factor  (1+c  )    by  the  Peetre 

inequality 

(l+n^)~^   <   2(l+c2)-l(i+|c-n|^) 

where  the  factor  (l+|c-ri|  )  can  be  compensated  at  the  expense 
of  differentiation  with  respect  to  k.  The  end  result  is  that 
P())  becomes  a  sum  of  terms  each  of  the  form 


gik-(c  n)[-^+j^2^-l~(^^^^^)^^^^^3^jj3i^ 

with  p  a  square-integrable  function  of  R-^xR^xr3  ^nd  ^ 
square-integrable  in  n .   This  shows  that  P(f)  e  3^  • 

We  now  state  the  second  theorem  of  Calderon  and 
Vaillancourt  [2]  in  the  following  simplified  version. 


THEOREM  ^.     The  operator  P,  defined  by  (17)  and  (I8)  with 
U  >  0,  can  be  extended  to  a  bounded  operator  f rom  5?f^  to;C 


13 


and  furthermore 

where  p   Is  the  constant  in  (l8)  and  the  constant  b  is 
independent  of  p(c,n,k). 

For  the  proof  of  this  theorem  we  refer  to  Ref .  2  .   It 
is  only  mentioned  that  the  extra  factor  (1+k  ) "^  in  (l8) 
enables  us  to  avoid  the  necessity  of  approximating  p(c,rijl<;) 
by  p .  of  compact  support  in  k  as  was  used  in  Ref.  2 . 

J 
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4.   The  Collision  Operator  (with  a  =  4) 

For  inverse-power  law  potentials,  the  scattering 
cross-section  has  the  simple  form 


|g|  I(|g|  ,0  )  =  lg|^  ^/^B(e,a)sin8 


and  for  the  Maxwell's  gas,  a  =  4  and  consequently 


(19) 


g|l(|g| ,6)  =  6(6,4)sine 


(20) 


which  is  independent  of  g. 


We  now  substitute 


.(c)  =  1  e^^*^$(k)d3k 


(21) 


into  (5)  and  obtain 


L(f)  = 


e^^'^p(c,k)$(k)d3k 


(22) 


p(c,k)  =  -exp(-tc^  -  ik'c}C[exp(^^+ik'c]  ]   (23) 


which  is  the  PsDO  form,  where  the  operator  C  is  given  in  (2) 
It  turns  out  that  some  of  the  integrations  in  (23)  can  be 
carried  out  explicitly  due  to  the  simplification  in  (20). 
The  result  is 


p(c,k) 
q(c,k) 


q(c,k)  +  r(c,k) 


^1T 


3(6,4)sin0[l-H(cose)]de 
3 


r(c,k)  =     B(e,4)sine  (2)  exp 
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f  1  2  1,  2  ., 
-p-c  -pk  -ik«c 


(24) 
(25) 

-  H(sin0) 
(26) 


de 


where   the    function  H  Is 


J. 

H(x)  =  r_i-— -l2  J 


{2-x^)- 


sln26 I cxk| 


L   (2-x^)   -i 


exp 


~x^(c^+k^+2ik«c)' 


(2-x2) 


(27) 


and  J  Is  the  zero-th  order  Bessel  function, 
o 

The  derivation  of  (24),  (25)  and  (26)  proceeds  as  follows 
We  write  (23)  as 


p(c,k,y)  =  expC-^(c+yk)^]C[expi(c+yk)^] 


(28) 


where  y  is  a  complex  variable  and  we  note  that  p(c,k,y=i)  =  p(c,k) 
But  first  we  consider  (28)  for  real  y.   For  real  y,  (28)  can  be 
integrated  but  by  using  the  method  of  Waldmann  [10]  with  his  X 
taken  to  be  -p-.   We  then  observe  that  for  fixed  c  and  k  the 
result  so  obtained  is  an  entire  function  of  y.   On  the  other 
hand,  for  fixed  c  and  k,  the  right-hand  side  of  (28)  is  also 
an  entire  function  of  y,  as  can  be  seen  from  (2).   Therefore 
we  can  use  analytic  continuation  and  set  y  =  i  in  the  resultant 
formulus  of  Waldmann  [10]  from  which  (24),  (25)  and  (26) 
immediately  follows. 

The  behavior  of  p(c,k)  can  be  studied  by  expressing  the 
dominant  contribution  near  8  =  ^tt  in  (25)  in  terms  of  confluent 
hypergeometric  functions,  which  we  shall  not  do  since  we 
shall  not  have  to  use  this  in  our  computations.   Suffice  it 
to  say  that  p(c,k)  is  C°°(R^xrJ)  and  satisfies  condition  (12) 


in  Sec.  3  with  a   =    ]x   = 
L  maps  ^  into  gf . 


1  _  1 
a    T" 


This  establishes  the  fact  the 
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As  we  mentioned  in  Sec.  3j  that  condition  (12)  with  a 
and  y  positive  makes  the  operator  (in  this  case,  N+L)  difficult 
to  treat.   The  difficulty,  arising  from  the  fact  that  diffe- 
rentiation with  respect  to  c  does  not  bring  down  the  growth 
in  c  in  (12),  is  caused  by  the  term  cxk  in  (25).   The  next 
section  shows  how  to  get  around  this  obstacle. 

Lastly,  we  note  that  the  singularity  of  3(G,a)  at  9  =  -^  it 
may  be  expressed  as 

6(0, a)  =  Bq(cos  e)"^"^^  +  BC9), 
|BCe)|  <  B^Cl  +  tan^  9)^, 
V  =  1/a . 
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5.   The  Associated  Operators 

The  key  to  getting  around  the  difficulty  raised  by  the 
term  cxk  mentioned  in  the  last  section  is  to  observe  that  if 
(f)  ( c )  is  a  wave  function  in  quantum  mechanics,  then  the  Fourier 
variable  k  corresponds  to  momentum  and  cxk  corresponds  to 
angular  momentum.   This  naturally  leads  to  the  consideration 
of  wave  function  with  fixed  total  angular  momentum;  i.e.,  we 
consider  (^    In  ^   having  the  form 


i(c)  =   Y^(c,z)iJ;(|c|  ) 


(29) 


and  denote  this  cj)  as  (^  c  J{l)>  where  Y  (c,z)  is  an  Jl-th  order 
spherical  harmonics  of  the  vector  c  with  z  as  the  polar  axis. 
The  Fourier  transform  $(k)  is  also  in^(£): 


$(k)  =  5(|k|)Yj^(k,z) 


(30) 


We    now   write    the    term   H(cos6  )exp  (ik 'C )    ir.    '25)    as 


H(cose)e^^*^    =    Q 


1+sln^e 


exp 


tco 
2(1+ 


s^e 


(1+sln   e) 


(c^+k^)] 


(31) 


Q  =  J, 


sin    9 

2 
"^l+sm    9 


IcxkllexpFfl-     ^°^    I    jlc-kl 
J  LI        l+sin    g''  -1 


1 

_-  CO 


y    a  J   ^u(C)P    (sin9)P    (cos®) 
^n   m  m+^  m  m '>         ^^  ■' 

m=0  -^ 


(32) 


where  P   are  the  Legendre  polynomials,  ^  = 2 | c | | k | sin9/ (1+sin  9), 
(h)  is  the  angle  between  c  and  k,  and  a   can  be  found  in 
Watson  [12  ]  . 
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It  is  noted  that  in  (32)  only  the  term  m  = £  can  act  on 
(30)  while  the  rest  of  the  terms  vanish  in  the  integration 
over  (H)  .   We  can  therefore  use  a  new  function  Q„ 


h  =   Pii(si'^9)c   ^a^J^+]^(^)Pn,(cos(H)) 


P  (sin6 )exp 


2sin  6 
1+sin^e 


ic  •!<; 


(33) 


where  Q.  has  the  same  result  as  Q  when  acting  on  c|)  of  the  form 
in  (30),  i.e.  on  (j)  £(/(£). 

Now  using  this  Q.  to  replace  Q  in  (25)  we  have  a  new 
symbol  q.(cjk;)  instead  of  q(c,k) 


q^(c,k)  = 


r^TT 


B(e,^)sinO[l-Hj^(cose)]d0 


(34) 


Hj^(x)  = 


3 


P„(y)    exp 


''-i-[(c2+k2)x2  +  2(l-y)2ic.k] 


(2-x2) 


K2-x")-i 
y  =  (1-  x^)^ 

With  a  similar  procedure  as  that  used  in  the  step  from 
(32)  to  (33)  J  we  can  construct  a  new  symbol  r.(Cjk)  in  place 


of  r (c ,k)  as 


rj^(c,k)  = 


^TT 


(35) 


6(0,4)siner(2)^exp(-|<;^-ik^-ik-c)  -  h^(sine)|de   (36: 


h„(sin  e)  =  H(sin  Q)        1'^  1  6  1  §'^5 


=  H^Csin  9)   0  £  e  <  ^TT, 

where  H^  is  given  in  (35),  and  H(sin  9)  is  given  in  (27).   We 
introduce 
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L,^    = 


e^^*''p^(c,k)$(k)d3k  ,  (37) 


p^(c,k)  =  q^(c,k)  +  r^(c,k)  .  (38) 

and  It  is  clear  from  the  construction  of  q,  and  r„  that 
L(p   =   L  (j)  If      (p   c  f^{i).      We  shall  call  L.  the  associated 

Jo  X^ 

operator  of  L.   If   cf)  ^  S^^'^  '      e.g.       ())  e  j{m)      with  m¥!i, 
then  Lcj)  ^   L„())  . 

The  main  advantage  in  dealing  with  the  associated 

operator  L.  is  that  the  derivatives  of  its  symbol  p  (c,k) 

behave  nicely  at  infinity.   Specifically,  we  have 

1     1 
Re  qj^(c,k)   >   a(c^+k^)^  +  ail   -  a^       (39) 

|rj^(c,k)|   <   a  (HO) 

_1 
|3'^3jjpj^(c,k)|   <   ad+k^)  ^,  I  n|  +  |m|  >  1    (^1) 


where  a  and  a,  are  positive  and  independent  of  i. 

We  first  go  to  (^0).   To  verify  (40),  we  split 
the  integral  in  (36)  into  two  parts,  one  for  0  <_  9  <_  jj-tt 
and  one  for  jf^   ±  Q   ±  jt^  •      In  the  interval  jt-tt  ^0  ^  p-ir ,  we 
can  factor  out  the  exponential  factor  exp  [-2<;  --^   -ik*c) 
which  ensures  the  boundedness  of  this  integral  for  all  c  and 
k;  in  fact,  it  ensures  the  decay  in  k  in  (4l)  for  the  Integral 
and  its  derivatives.   In  the  other  interval  0  <_  6  ^  jt-tt  , 
6(9,4)  is  finite  and  (40)  follows  from  this. 


20 


To  prove  (4l),  we  only  note  that  q  (x,k)  and  h  (sin  Q) 


I 


I 


2    2 
only  contain  c  ,  k   and  c-k  and  therefore  differentiation  with 

respect  to  c  and  k  can  be  easily  carried  out  and  (4l)  can  be 

established  almost  straightforwardly. 

To  prove  (39),  we  split  (3^)  into  three  terms 


u 


q^(c,k)  =  Q-^  +  Q2  +  Q3, 

6(e,4)sin  e  [1  -  Pj^(y)]e"'°  de, 

-r 


-Ft 


J  0 

1. 


«2  = 


%- 


6(e,  4)sin  Q[l  -  e    ]d9. 


2TT 


6(G,'^)sin  G  Pj^(y)e   °[l-r^]dQ, 


(42) 

(43) 

(44) 

(45) 

where  the  functions  F   and  T,  are 

o      1 


r^  =  I  x2(c2+k2)/(2-x2). 


r^  =  [2/(2-x^)]3/2exp[-i(l-y)^c-k/(2-x2)] 


(46) 
(47) 


We  first  note  that  Q^  is  bounded  in  all  c,  k  and  I. 

For  Q„,  we  can  introduce  a  new  integration  variable 

,    2^,  2sl/2 
t  =  (c  +k  )    X 


(48) 


and  show  that 


Re  Q2  >  2a(c^+k^)^^^  "  |  ^i 


(49) 


for  some  positive  a  and  a,  . 


On  the  other  hand,  for  Q.,  we  can  use  the  formula  of  Hilb; 


see  Szego  [9]  : 


Pj^(sin  G)  =  Jq[(£+|)(^-G)]  +  R^(G)  , 
R^(0)|  <_   a|G|  . 


(50) 
(51) 
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Now  we  substitute  (50)  into  (43),  Introduce  a  new  in- 
tegration variable 

t  =  {1+   |)(|^-G), 

and  use  the  inequality 

r^  <  ^  t2(c2+k2)/(£+  h^. 

o  —  2  2 


As  a  result,  we  obtain  for  Q-, 


Re  Q^  >_   6^£ 


t  ^-2^[1-J^(t)]exp[-  |(c2+k^)/(£+  |)^]dt 


+  b(c,k) 

where  b(c,k)  is  uniformly  bounded  for  all  c,k  and  £.  In  the 

above  estimate  we  have  used  the  explicit  form  of  the  singularity 

of  6(0,x)  at  G  =  p-TT  given  at  the  end  of  Sec.  4.   An  alternative 


way  of  writing  the  last  estimate  is 


Re  Q-,  ^  al 


2   2     2 
a^,   c  +k  <  Z 


(52) 


2   2      2  1/2 

while  for  c  +k   >  i      we  can  use  (49)  to  get  a  i  factor. 

In  any  case,  (49  )  and  (52)  together  with  the  boundedness  of  Q- 
lead  to  the  desired  result  (39)- 
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6.   The  Eigenvalue  Problem 


We  now  consider  the  eigenvalue  problem 


A({)  =  A^<{.  =  (N  +L^)(}) 


(53) 


where  L.  is  given  in  (37)  and  N  is  a  large  but  fixed  positive 
number.   What  we  like  to  do  is  to  invert  the  operator  Ag.      In 
doing  this,  we  essentially  follow  the  general  procedure 
described  in  Kumano-go  [  7],  making  use  of  the  four  theorems 
in  our  Sec.  3. 

First  we  rewrite  the  operator  A„  as  a  PsDO  with  multiple 
symbol  in  the  form 


\^   = 


ik*  (c-n)   /  ,  \  a  /  ^^3  j3, 
e       "a^(c,k)<})(Ti)d-^nd^k  , 


(5^) 


aj^(c,k)  =  N  +  Pj^(c,k) 


(55) 


The  symbol  a.(c,k)  may  be  expressed  as 


a^(G,k)  =  a^(n,k)  +  b(c,n,k) 

1 
b(c,n,k)  =  (c-n)'    [3^Pjj^(ri+x(c-ri)  ,k)  ]dx 

o 


(56) 
(57) 


Now  substituting  (56)  and  (57)  into  (5^)  and  integrating 
by  parts  yield 


A^<p   =   A^<t>   +   T()) 


(58) 


R  R 

where  A„  is  the  reverse  operator  [  7  ]  of  A  and  A   and  T  are 

given  by 
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aJ*  = 


f  ei^-(^-n)a^(n,k)((.(n)dndk 


(59) 


T<t.  =  I  I  e^^*^^-^^(c,n,k)<j,(n)d3nd3k 


(60) 


t(c,ri,k)  =  i 


[8j^-8^aj^(n+x(c-n)  ,k)]dx 


(61) 


It  is  noted  that  In  the  Integration  by  parts  we  have  to 
follow  the  procedure  of  Kumano-go  (see  Ref.  7  ,  pp ■  ^21-422), 
using  the  fact  the  (fi  e  ^.   According  to  (4l),  the  multiple 
symbol  t(Cjri,k)  satisfies  the  condition  in  Theorem  4  and  conse- 
quently T,  which  maps  jjf  into  ad  ,  can  be  extended  to  a  bounded 
operator  from^^  to;56^  with  li  Til  independent  of  N. 

Next  we  introduce  the  operator  B  defined  by 


B(J) 


=  I  e^^-^[a^(c,k)]-l$(k)d3k 


(62) 


N 


where  Re  a„(c,k)  >  -,  if  N  is  large  enough,  according  to  (39) 
^        8 

and  (40).   Now  Theorems  1  and  2  imply  that  B   is  a  compact 

2 
operator  on ^     and 


Bll   <   b/N 


(63) 


for  some  constant  b 


We  then  consider  the  composite  operator 


BA^  =  BA^(})  +  BT(j) 


(64) 


and  write  the  operator  BA.  as 


BA^cj)  = 


gik.(c  'n\si^(T),k)/a^{c,k):i<p{^)d^nd\ 


(65) 
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Now  substituting  (55)  for  a,(n,l<:)  in  (65),  we  obtai 


n 


BaJ(()   =   ())(c)    +   Y(})  (66) 

YO)  =  I  I  e^^^^-^^y(c,n,k)<{,(n)d3nd3k  (6?) 

1 
y(c,n,k)    =   i   J    dx9j^-[a^^(c,k)3^a^(n+x(c-n),k)]  (68) 

o 

From    (^1),    we    see    that    y(cjrijk)    also   satisfies    the 
condition   in  Theorem   4   and 

IIYII       <      b/N  (69) 

for  some  constant  b,  independent  of  £. 

Now  summarizing  the  result,  we  have 

BA^*  =  (1  +  Y  +  BT)4)  (70) 

and   the   bounds    in    (63)    and    (69)    show   that 

II Y   +   BTII       <      b/N  (71) 

for  some  b,  from  which  it  is  concluded  that  if  N  is  large 
enough  then  (BA  )    exists  and  is  bounded  in ^  .   Therefore 
the  eigenvalue  problem  (53)  may  be  written  as 

^  =  (BA^)-^B<},   =   A~^(^  (72) 

2      2  -1 

and  since  B  is  compact  from,;^^  to  jz^  and  (BA  )    is  bounded 

—1  2      2 

it  follows  that  A.   is  a  compact  operator  from^f,  to  si  . 
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Next,  we  wish  to  show  that  as  il  ^  oo,  we  have  ||A,  ||->  0. 
In  this  connection,  we  first  note  that  because  of  (39) 'we  have 

|[a^(c,k)]"^|  -  0 

as  £  ^  00.   This,  by  Theorem  1  of  Sec.  3,    shows  that  ||b||->  0  as 
£  -»■  oo.   On  the  other  hand,  the  constant  b  in  (71)  can  be  made 
Independent  of  £,  and  therefore 


A^^ll  ^  0  (72) 


as  £  -»•  00. 


-1  2     2 

The  fact  that  A.   is  a  compact  operator  from  'ji     to  ^ 

together  with  the  result  in  (72)  leads  to  the  following 

conclusions. 

(1)  A    is  a  compact  operator  from  X   to  ^  . 

(2)  The  eigenfunction  of  A  are  all  ,^  -functions.   Since 

p 
A  =  N+L,  the  eigenfunction  of  L  are  all  ^  -functions. 

(3)  The  spectrum  of  L  is  purely  discrete. 

2  2 

(4)  The  %   -eigenfunctions  of  L  are  complete  in  ^  ,  since 

L  is  symmetric  and  non-negative. 

(5)  Further  computation  shows  that  the  eigenfunctions  of 
L  in  fact  belong  X,o  A   where  j^  is  the  set  of  C  -functions  which, 
together  with  their  derivatives,  decays  at  infinity  faster  than 
any  inverse  powers  of  c .   This  will  be  proved  in  Part  II  for 
the  general  case. 

(6)An  alternative  way  of  stating  the  first  conclusion  is 

^1  0  0 

that  L   is  a  compact  operator  from  ^  to  X  • 
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7.  Remarks 

Our  analysis  relies  on  two  main  tools.   One  Is  the  theory 

of  pseudo-differential  operators  and  the  other  the  use  of  the 

associated  operator  A..   In  constructing  A.,  the  major  step 

Is  the  step  from  (32)  to  (33).   Since  only  the  term  with  m=  I 

In  (32)  need  be  kept  unchanged,  there  Is  considerable  freedom 

In  altering  the  other  terms  with  m  7^£.   One  way  would  be  to 

set  all  of  them  to  zero,  leading  to  the  result 
_1 
Q^  =  C  ^A^J^^3^(OP^(slne)Pj^(cos  ©)  . 

If  we  do  this,  we  would  end  up  with  a  one-dlmenslonal 
equation  with  |c|  as  the  Independent  variable,  which  Is  what 
Wang  Chang  and  Uhlenbeck  have  In  Ref.  11.   Then  we  either  have 
to  construct  the  elgenfunctlons  or  have  to  deal  with  (Instead 
of  three-dimensional  Fourier  transforms)  one-dlmenslonal  Hankel 
transforms  for  which  there  exists  no  pseudo-differential 
operator  theory.   This  consideration  leads  us  to  the  construc- 
tion in  (33)  which  retains  the  pseudo-differential  operator 
formalism. 
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Part    II 


General  Repulsive    Inverse-Power   Potentials 
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1.   Introduction 

This  paper  studies  the  eigenvalue  problem  of  the  collision 
operator  in  the  linearized  Boltzmann  equation  with  a  repulsive 
inverse-power  intermolecular  potential  of  the  form 

V(r)  =  ar"",  a  >  2  (1) 

where  r  is  the  intermolecular  distance. 

The  purpose  of  this  paper  is  to  show  that  the  collision 

operator  L  [see  eqn.  (12)]  has  a  purely  discrete  spectrum  and 

v^2 
its  eigenfunctions  exist,  are  X   -functions  and  form  a  complete 

2 

orthonormal  set  in  ^   .      The  result  is  obtained  by  showing  that 

L~   is  a  compact  operator  from  ^  to  %.,   where  ^j_  is  the  sub- 

2 
space  of  ^  orthogonal  to  the  five  eigenfunctions  of  L  for  the 

2 
eigenvalue  zero,  and  ^  are  the  square  integrable  functions. 

The  method  we  use  here  follows  closely  the  procedure 
carried  out  in  the  first  part  of  this  study  [4]  (hereafter 
referred  to  as  Part  I),  for  the  special  case  of  a  =  4,  i.e. 
the  Maxwell's  gas.   It  is  suggested  that  Part  I  be  read  as  a 
preliminary  since  the  present  report  makes  frequent  use  of  the 
results  given  in  Part  I. 

In  the  next  section  we  consider  the  operator  A  =  N+L, 
where  N  is  a  large  positive  number  and  L  is  the  symmetrized 
collision  operator  in  (12),  and  we  construct  the  associated 
operator  A,  in  Sec.  3.   We  use  the  formalism  of  pseudo- 
differential  operators,  [see  Sec.  3j  Part  I].   The  properties 
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of  the  symbol  of  the  operator  A„  is  investigated  in  Sec.  4. 
These  properties  (in  terms  of  bounds  for  the  symbol  of'  A. 
and  its  derivatives)  are  then  used  in  Sec.  5  to  show  that 
A~  exists  and  is  a  compact  operator  from^c  to ^  .   From 
this  we  conclude  that  the  spectrum  of  L  is  purely  discrete 
and  its  eigenfunctions  are^^  -functions  and  they  are  complete 
in  ^  .   Further  computations  show  that  the  eigenfunctions  of 
L  are  not  only  ^  but  also  belong  to  jjf ,  the  rapidly  decreasing 
C  -functions.   These  computations  are  sketched  in  Sec.  6, 
followed  by  some  general  remarks  in  Sec.  7. 
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2.   Collision  Operators 

Consider  the  collision  operator  C  in  the  linearized 
Boltzmann  equation 


C[f]  = 


K'(c,n)[f(n)-f(c)]exp(-n^)d3n 


+  I  [K^(c,n)-K^(c,n)]f(n)exp(-n^)d3n   (2) 


where  kernels  K' ,  K'  and  K-,  correspond,  respectively,  to 
terms  involving  what  is  conventionally  denoted  as  f ,  f'  and 
f-i  in  the  collision  integral.   These  kernels  may  be  written 


as  [1] 


K'(c,ri)  =    sine  G(tane  )sec^eo  (  |  g  |  secB  ,  6  )d8 


K^(c,n)  = 


f  ?    ,  , 

sine  G(cote)csc-^ea(  I  g|  csce,e)de 


(3) 


(4) 


K^(c,n)  = 


sine  a( |g| ,e)de 


(5) 


=  n  -  c 


G(y)  =  I^(2|cxg|y)exp(-y  g  ) 


(6) 
(7) 


where  we  have  used  c,  n ,  g  to  denote  vectors  and  |g|,  |c|, 
I n I  to  denote  scalars  and  the  e-integration  in  (3)-(5)  is  from 
0  to  J^TT .   The  function  a(|g|,6)  is  related  to  the  differential 
cross-section  and  for  inverse-power  potentials  given  in  (1), 


o(|g|  ,0)  =  |gr6(e,a) 


(8) 


where  y=  1  -  ^^ ,    v= 1/a ,  y    >    -1      for   a  >  2 ,  and  3 ( 6 , a ) 
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diverges  as  (^tt-O)    ^  as  0  ->  ^tt  .   This  property  of  6(0, a) 
may  be  expressed  as 

e(0,a)  =  3Q(cos0)^"'*^(tane)^"^^+  B(tan0)  (9) 


|B(tan0)|   <   b(l  +  tan^0)^ 


(10) 


which  will  be  found  useful  later.   We  may  just  set  ^q  =  1  on 
account  of  the  linearity  of  the  eigenvalue  problem.   Throughout 
this  paper  we  shall  use  b  to  denote  positive  constants,  not 
all  of  them  necessarily  the  same. 

It  is  expedient  to  use 


())(c)  =  f(c)exp(-^c^) 


(11) 


and  consider  the  operator 


L(l)  =  -exp(-^c^)C[4)(c)exp(^c^)]  . 


(12) 


which  can  be  shown  to  map  jjO  into  ^     (this  was  discussed  in 
sec.  3  of  Part  I)  and  is  symmetric  and  non-negative.   Conse- 
quently L  can  be  considered  as  a  non-negative  self-adjoint 
(but  unbounded)  operator  in  ^  .   Our  main  task  is  the  study 
of  the  eigenvalue  problem 


Xcf)  =  (N+L)(})  =  A()) 


(13) 


where  N  is  a  large  but  fixed  positive  number. 

To  write  A  in  a  pseudo-differential  operator  (PsDO) 
form,  it  suffices  to  use  the  Fourier  transform 


<})(k)  =  (2Tr) 


-3 


-ik  'c  ,  /  \  ,3 
e      ({)  (c  )d  c 


d'^) 
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in  terms  of  which  Lij)  becomes 


Lcf)  =  I  e^^*^p(c,k)$(k)d3k 


(15) 


p(c,k)  =  q(c,k)  +  r(c,k)  +  u(c) 


(16) 


q(c,k)  =  I  d3gK'(c,n)[l-e^^*S]exp(-^c2-^n^)  (17) 


r(c,k)  = 


d3g[K^(c,n)  -  K|(c,n)]exp(-^c2-]^n^+lk-g)   (18) 


u(c)  =  j  d^gK-  (c,n)Cexp(-^n^)-exp(-]^c^)]exp(-3^n^)  .   (19) 

Equation  (15)  gives  L(j)  in  the  form  of  a  PsDO  and  p(c,k) 
is  the  symbol  of  the  PsDO  L.   Using  the  analyses  in  Sec.  3,  we 
can  show  that  p(c,k)  satisfies  condition  (12)  in  Part  I  for 
a  =  n  =  V,  and  consequently  L  maps  /^  into;*^,  according  to 
Theorem  1,  Sec.  3,  of  Part  I.   But  we  shall  not  carry  out  the 
detailed  estimates  since  we  shall  not  have  to  use  the  estimates 
of  p(c,k)  and  its  derivatives  in  our  subsequent  calculations. 

The  operator  A  then  becomes 


AcJ) 


'  e^'^*^a(c,k)$(k)d^k  , 


(20) 


a(c,k)  =  N  +  p(c,k)  . 


(21) 
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3.   The  Associated  Operators 

As  it  stands,  the  operator  A  in  (19)  is  difficult  to 
invert.   This  difficulty  has  been  discussed  in  the  special 
case  of  a  = 4  (the  case  of  Maxwell  molecules)  and  it  can  only 
become  worse  for  the  general  case  of  a  ^  >4 .      The  way  out  of 
this  is,  as  we  have  done  in  the  case  of  Maxwell  molecules,  to 
construct  the  associated  operator  A.  (see  Sec.  5,  Part  I). 
The  construction  for  such  an  A„  in  the  case  of  general 
inverse-power  potential  with  a^H,    unlike  the  case  of  Maxwell 
molecules,  is  somewhat  complicated  and  a  brief  explanation 
would  be  quite  helpful.   This  may  be  given  in  the  following 
heuristic  arguments. 

The  main  difficulty  is  in  the  symbol  q(c,k)  which 
arises  from  the  term 

C^Ef]  =  J  d3nexp(-n^)  I  dcf  J  de|g|^6(0,a)[f  (c'  )-f(c)] 
in  C[f].   However,  we  observe  that  for  large  c 

C^[f]  'x.  |c|^  I  d3nexp(-n^)  J  dcj)  [  dOg  (  6  ,a)  [f  (  c  '  )-f  (c  )]   (22) 

which  has  the  same  form  as  the  collision  operator  for  the 
Maxwell  molecules  except  for  the  factor  B(0,a)  instead  of 
3(6,4)  and  the  extra  |c|^,  both  of  which  can  be  easily 
incorporated  into  the  treatment  in  Part  I.   It  is  therefore 
natural  that  we  shall  explore  this  large-c  behavior.   This 
will  now  be  done  in  a  rigorous  manner  in  the  following. 

We  shall  now  construct  the  associated  operator  A. 
in  four  steps. 
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First    we    rewrite    q(c,k)    as 


q(c,k)    =    s(c,k)    +   v(c,k) 


v(c,k)    =    q(c,k)    -    s(c,k) 


kir 


S(c,n)    = 


3, 


slneG(tane)3(e,a)sec-^0(sln0) 


,-l  +  i|v 


de 


(23) 
(24) 


s(c,k)    =    |c|Y   I   d3gS(c,n)[l-e^^'S]exp(-i^c2-]^n2)  (25) 


(26) 


where  G(tan0)  Is  given  in  (7)  and  we  note  that  the  kernel 

S(c,ri)  is  precisely  the  kernel  for  (22)  except  for  the  extra 

4v-l 
factor  (sine)     .   Since  the  singularity  of  3(9, a)  is  at 

6  =^TT,  the  extra  factor,  being  equal  to  1  at  6  =  ^tt  ,  does  not 

have  any  particular  significance  and  is  added  here  merely  for 

later  computational  convenience. 

The  next  step  is  to  introduce  a  C  -function  Q{c)   with 
the  property  that  for  some  c   >  1,   c(c)  =  1   for  |c|  <  c  , 


^(c)  =  0  for  |c|  >  2c  ,  and  0  <  c(c)  <_  1  for  c   <  |c| 
and  we  write 


<  2c  , 
o' 


q(c,k)  =  c(c)q(c,k)  +  [1-^ ( c ) ] [s ( c ,k )+v( c ,k ) ] 


(27) 


So  far,  nothing  has  changed  in  q(c,k);  we  merely 
rewrite  q(c,k)  in  a  form  from  which  it  is  convenient  to 
construct  the  new  operator  A„. 

The  third  step  is  to  rewrite  s(c,k)  in  (25)  as 


s(c,k)  =  I c 


d^n 


^j.jno»co\r  ik'c   ik'c'-i 
d(|)d63*  (e)  [e    -e     ] 

X  exp[-  |c^  -  n^  -  ik-c  +  |(c')^]  (2^ 
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which  is  Just  the  type  of  integral  that  appears  in  the 
collision  integral  for  the  Maxwell  molecules.   Here  3*  (6)  is 

3*(0)  =  6(0,a)(sine)^^"^  .  (29) 

The  integrals  in  (28)  can  be  evaluated  by  using  the  formulus 
of  Waldmann  [5]  in  conjunction  with  analytic  continuation  as 
we  have  shown  in  Sec.  5,  Part  1.   The  result  becomes 


s (c,k)    =    I c 


B*(e)(2t) 


3/2 


X    ex 


p(-3^tcos^ec^)ll  -  JQ(tsin2elcxk|  ) 
X    exp{-^tcos^e(k^+2ic'l^}Jde 


(30) 


t  =  (1  +  sin^e)"-^ 


(31) 


We  can  now  repeat  the  procedure  in  Sec.  5  of  Part  I 
to  construct  a  new  symbol.   In  fact,    we  consider 


r 


e^^*^s(c,k)$(k)d3K 


(32) 


for  (j)  e  ^  { I)  ,    viz  . 


(})(c)  =  Y^(c,z)4;(|c 


(33: 


where  Y.(c,z)  is  the  £-th  order  spherical  harmonics  for  c  with 
polar  axis  z.   Then  $(k)  is  also  in  ^(  £ ) ,  and  we  follow  the 
same  step  of  Part  I  [see  eqns.  (32)  and  (33)  of  Part  I]  and 
obtain  a  new  symbol 


4o 


s^(c,k)  =  |c 


3^71- 


B(e,a)(slne)^^(2t)3^^exp(-^tcos^ec^) 


X  [1  -  P^(slne) 'expC-r)] 


(3^^) 


r  =  l^tCk^cos^e  +  2(l-slne)^lc*k] 


(35) 


where  t  =  p —  .   We  note  that  for  <t>    z  ^ {I)    in  (33), 


(1+sin  e) 


we  have 


f  .lk-c_,.  ,.^^,,.^.3v  =  \    ^il<*c,  (^   ^.\l(^^r^%.  (36) 


e^^  ^s(c,k)(})(k)d^k  =  j  e^^  ^Sj^(c,k)4.(k)d^k  . 

For  a  detailed  discussion  of  this  point  we  refer  to 
Sec.  5  of  Part  I. 

Now  the  final  step  is  to  construct  a  new  symbol 
q^(c,k)  from  (27)  as 

q^(c,k)  =  eC(c)Sj^(c,k)  +  [  1-C  (c  )  ]  [s^  (c  ,k)  +  v(c,k)] 

(37) 
+  C(c)[q(c,k)  -  es  (c,k)] 

where  we  have  changed  two  s(c,k)  terms  to  Sjj^(c,k)  with  others 

remaining  unchanged  as  in  (27). 

With  this  new  symbol  q  (c,k),  two  new  operators  L 
and  A   may  be  introduced  as 


p  (c,k)  =  q  (c,k)  +  r(c,k)  +  u(c) 

A/  A/ 

L  (j)  =    e"^  *"^p  (c,k)$(k)dk  , 


A. 4) 


'  e^^*^aj^(c,k)$(k)dk  , 


a^(c,k)  =  N  +  Pj^(c,k) 


(38) 

(39) 
(40) 
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where  the  main  point  to  note  is  that  A<i>   =   AA   for  <p   €  ^(£) 
in  (33).   It  should  be  pointed  out  that  since  L  is  isotropic 
all  eigenf unctions  of  A  are  inj^(A)  for  some  £. 

What  we  shall  show  is  that  A~   exists  and  is  a 
compact  operator  from ^  to  SC  .   This  we  shall  prove  by 
using  the  theorems  of  PsDO  given  in  Sec.  3,  Part  I.   For 
this  purpose  some  properties  of  the  symbols  q^(c,k),  r(c,k) 
and  u(c)  have  to  be  established;  this  is  undertaken  in  the 
next  section.   These  estimates  can  also  be  used  to  show 
that  II A^  II  ^-  0  as  «,  ^-  oo,  which  leads  to  the  fact  that  A""^ 
is  compact  from  X  to  ^  . 
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H.      Properties  of  the  Symbols 

The  purpose  of  this  section  Is  to  establish  some  bounds 
for  the  various  symbols  and  their  derivatives  with  respect  to 
c  and  k.   These  estimates  Involve  considerable  computations. 
But  they  are  of  critical  Importance  to  our  theory  and  some  of 
them  are  not  at  all  very  obvious.   For  these  reasons  we  shall 
present  here  the  major  steps  In  the  proofs  of  these  estimates. 
It  is  recalled  that  throughout  this  paper  we  shall  use  b  to 
denote  positive  constants,  not  necessarily  all  of  them  the 
same  constant. 

Let  us  first  summarize  the  main  results  of  this 
section : 

N  +  Rea„(c,k)   >   ^N  +  b (1+c^ )^^( c^+k^ )^  +  bi^^  (42) 

|3^a^(c,k)|   <   bd+k^)^  ,     lm|  >  1  (43) 

|a^9jja^(c,k)|   <   b(l+k^)   ^  ,   |n|  >  1  and  |m|  >0   (44) 
where  the  constants  b  in  (42)  -  (44)  are  independent  of  £ 

and    V,  =  v(l-2v)/(l+2v)  >  0  and  Re  denotes  the  real  part 

and  we  have  used  the  abbreviation 

m,  m„  m-j 
9"^  =  9  ^3  ^9  ^  ^  (45) 

|m|  =  m-,  +  mp  +  m^  ,  (46) 


wh 


ere  the  c,  1  =  1,2,3  are  the  three  components  of  the  vector  c, 


Similar  notation  is  also  used  for  8,  . 

Theresults  in  (42),  (43)  and  (44)  are  all  that  are 
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needed  for  the  analyses  in  the  next  section.   We  now  embark 
on  the  proof  for  these  estimates  by  considering  the  terms 
u(c),  r(c,k),  c(c)q(c,k),  v(c,k),  Sj^(c,k)  and  ^(c)s(c,k)  in 
(17),  (18),  (19),  (24),  (3^)    and  (30),  one  at  a  time.   Those  not 
interested  in  the  details  of  the  proofs  can  skip  the  rest 
of  this  section. 

(A)  Estimates  for  u(c) 

To  obtain  the  desired  estimates  for  u(c),  we  shall 
reduce  the  five-fold  integral  in  (19)  into  a  one-fold  Integral 
of  confluent  hypergeometric  functions.   To  this  end,  we  denote 
i)   as  the  angle  between  c  and  g  and  consider 

(  2  2    2 

U  =    sinij^d^j*!  (2  |  cxg  1  tan  e)  [exp  (-n  )-exp(-Jin  -ijc  )] 

which  can  be  Integrated  out  in  terms  of  Bessel  functions  of 

,    1 
order  ^   as 

U  =  (2^)^[I^(Ci)c-^  -  Ij^(^2)^2^^  ^^^) 

^1  =  2|c| |g|/cose  (48) 

^2  =  2|c|  lg|  (tan^e  +^)^  (49) 

by  using  the  well-known  addition  theorem  of  Bessel  functions, 
see  eqn.  (9)  on  page  370  in  Watson  [6], 

Now  substituting  (47)  into  (19)  and  using  the  table 
for  Hankel  transforms  [3]  yields 

u(c)  =  b    de3(e,a)(sece)^  ^sinG  Wexp(-c'^)  (50) 
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where  b  is  a  positive  constant  whose  exact  value  is  of  no 
importance  to  our  analysis,  and.  W  is 


W=  {cosQ)^'^'^^F^[^+^y,^,  +c^]    - 


■3  ^  1„  3 


niQ) 


,  ]^  _  2  -  ^2Y 

=  [tan^e  +  ij^ftan^e  +  i) 


(51) 


(52) 


-1 


6  =  i  [tan^  e  +  i-}    , 


(53) 


where   F   is  the  confluent  hypergeometric  function. 

To  obtain  estimates  for  u(c)  in  (50)  and  its  deriva- 
tives, it  is  convenient  to  regroup  W  in  (51)  as 


W  =  W-^  +  W2 


W-L  =  [(cose)3  +  Y-n(6)]^F^(|+^Y,|>c^] 

W2  =  f^(e)^^F^||+|Y,|,c2J  -^F^^|  +  |y,|,c2-6c2J 


(54) 

(55) 
(56) 


Substituting  (5^)  into  (50)  and  integrating  by  parts 
with  0  for  the  Wp  term  yields 


u(c)  = 


u^(c)  +  U2(c) 


u-j^(c) 


"  ^o  2Pi[^  +  ^y>i>^^)^^p(-^^^ 


(57) 
(58) 


U2(c) 


=  be' 


^^^   d6  ~ 


de 


d6 


3(9)^F-l(|+  ^Y,|,c^-6c^)exp(-c^)    (59) 
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where  a   is  the  Integral  in  (50)  of  the  bracketed  term  in  (55) 
which  may  be  negative,  and  6(6)  is 


e 

B(e)  =  (I  +  y/3)  I  d0sine(sece)3'^^3(0,a)f2(e)  .      (60) 

o 

The  behavior  of  B(9)  can  be  obtained  from  (52),  (9) 

and  (10)  and  may  be  expressed  as 


3(9)  = 


-1   r-v  ,  ^rr-v+^-i 


l+2v 


5-v  +  0(6-^"^) 


(61) 


where  6  is  given  in  (53). 

For  large  |c|,  the  confluent  hypergeometric  function 
has  the  asymptotic  behavior  [2] 


lF-L(|+i-Y,|,c2-6c2) 


2^  "^ 


r(|) 


Lr(i^^Y)J 


[c^(l-6)]^exp(-6c^+c2)   (62) 


Substituting  this  into  (59),  we  have  for  large  c 


,    s    "^       h        11  Y+2v 


r(|)  -^ 


r  (|4y) 


e   X  dx 


(63) 


and  on  the  other  hand,  from  (58)  we  have  for  large  c 


u^(c)l   ^   b|c|^ 


(64) 


With  the  result  in  (63)  and  (64),  the  following 
estimates  are  obtained 


u(c)   >   b 1 c 


Y+2v   1 


N 


u(c)l   <  bd+c^)^^"*"^ 


(65) 
(66) 
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for  all  c  if  N  is  taken  large  enough. 


The  estimates  for  the  derivatives  of  u,  and  Up  can 
easily  be  obtained  as  follows.   First  we  use  [6] 


^F^(|  +  |Y,|,+c^)exp(-c^)  =  iF^(-|y,|,-c^) 


(67) 


9z  1^ 


1„  5 


F^(-i-Y,f,z}  =  ^F^(l-^Y,|>z)(-f) 


(68) 


1,.  5 


1 


Fl(l-^Y,|,z)  I  1  b|z 


^^ 


as   z  ->  -00  , 


(69) 


which  shows  that  each  differentiation  with  respect  to  c  brings 
down  a  power  of  |c|  as  c  ^  oo.   Similar  argument  applies  to 
Up(c),  and  we  therefore  have 

>   1  (70) 


3%(c)|   <   b 
c     ' 


m 


The  estimates  in  (65),  (66)  and  (70)  are  the  desired 


results 


(B)  Estimates  for  r(c,k) 

To  obtain  estimates  for  r(Cjk)  in  (I8)  and  its 
derivatives  with  respect  to  c  and  k,  we  first  split  the 
e-integration  into  two  ranges,  0  <_  6  <_  9  ,  and  9   <  9  <_  p-TT, 


<  j-n ,    such  that  |cot9|  <  jr-  in  9   <_ 


<_  p-TT .   We  then  write 


r(c,k)  =  r-|^(c,k)  +  r2(c,k)  +  r2(c,k) 


r-|^(c,k) 


d^g    d9sin9B(9,a)|g|^Cl-(csc9)^"^^G(cot9)] 


X  exp  [-pC 


12   12,. 


+  ik'g) 


(71) 
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u 
r2(c,k)  =  J  d^g  J  de  6(e,a)(csce  )^"^^  G(cot6  )|g|^ 


X  exp  (-  2  c  -  2  n  +  ik'g) 


r^(c,k)  =    d-^g    de   sine  3(6, a) 

0 


^12    12,.,    V  I  lY 
X  exp(-  p  c   -  p  n  +  ik'g)ig|' 


In  (73)  we  note  that  the  function   6(6, a)  sine   Is 
absolutely  Integrable  In   0  1  6  <_  e„  ,  and  consequently 

|r2(c,k)  I  £  b  exp(-  2"  c  )    dn  exp(-n  )  |  c-n  |^ 
<  b(l+c^)^/^  exp(-  I  c^) 


(72) 


(73) 


(74) 


On  the  other  hand,  we  can  write 

l|k|g 


Ik'g    1 
1  e    °  =  


[^ 


-    1) 


(75) 


|k|  8g„ 

and  Integrating  by  parts  with  respect  to  g   ,  the  component 
of  g  parallel  to  k.   We  then  obtain 


r^(c,k)|  <  b(l+c^)^'^^(l+k^)  ^^^  exp(-  ^  c^) 


3 


The  derivatives  of  r^(c,k)  with  respect  to  c  and  k 


(76) 


can  be  estimated  similarly  and  the  result  may  be  written  as 

(77) 


ti'"  ^n  ^  .   ,  N  I  ^  v,/-,  ,,  2s-l/2,,  ,  2>-l/2 
d   9,  r^(c,k)|  <_   b(l+k  )     (l+c  ) 


For  r-,(c,k)  and  r„(c,k),  we  have  to  use 


>m  ^n 


2  X  m  +  n 


^c  ^g  ^o^^lcxglcot  e)|  <  b(l+c^)''"l^l"l  lQ(2|cxg|cot  e)   (78) 


obtained  from  the  series  for  1„  by  termwlse  differentiation, 


-4  8- 


Now  for  r-,  (c,k),  we  note  that  in  the  range 


1 


<  6  <_  i-TT ,  I  cote  I  _<  ttj  and  consequently 


9'^9'^G(cot9 
eg' 


2sM-, 


...|1- 


I   <   b(l+c^)^"^lj2|cXg|^) 
<  b(l+c^)^exp[|lc|  In  I)  . 


(79) 


We  first  Integrate  (71)  by  parts  with  respect  to  9, 


usxng 


e 


3(0) 


de  sine3(e,a)  , 


e 


we  obtain  the  estimate 


r^(c,k)|   <   b 


d3g^(l+g2) 


^7T 


3(e) 


(80) 


X  exp(-^  -  2-'l  +  2"!  "^  I  '  "^  I -^^^ 

1+^ 
<_     b[l  +  c  ]    exp  (-• jT-c  ) 


(81) 


Now,  using  the  same  procedure  of  partially  integrating 
with  respect  to  g,,  as  in  the  estimate  for  rT(c,k),  we  have 

2+1 

|r^(c,k)|   <   b[l+c^]   ^exp(-^c^)  •  (l+k^)~^ 


l+k2)-^(l+c^)"^ 


(82) 


<  bd+k'^)  "2(1+0  ) 


And  the  derivatives  of  r^(c,k)  can  be  estimated  in 
the  same  fashion  by  using  (79)  and  the  result  is 


a'^aj^r,  (c,k)l  <  b(i+k2)-^(i+c2)-^ 

c  k  1   '   '   — 


(83) 


for  Iml  >  0  and  In   >  0. 
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Finally  we  turn  to  rp(c,k)  which  is  the  most  difficult 
part  of  r(c,k).   We  first  introduce  x  =  |g|cot9  and  note  that 


3(e,a)(csce)^^^de  =  B(x,|g|)dx  , 

[l+(x/|g|)^] 
|B(x,|g|)|   <   — 


(84) 
(85) 


which  is  readily  seen  from  the  fact  that  B(9,a)  is  bounded 


for  0  <  6  <  0   with  9  <   rrT\ , 
—   —  o       o        d 


We  now  write  (72)  as 


|r2(c,k)|  < 


d^g 


-1 


dx|gr^[|gr  +xM 


12      12 


X   i^(2|    cxsin(H)  |)exp(-2-c    -  ^-n    -X   ]    ,      (86) 


where  (B)  is  the  angle  between  c  and  g. 

It  is  useful  to  make  the  following  observations: 
(1)   The  contribution  to  (86)  from  the  part 

-1  "I   O 

|sln(H)|    <      jr  is   bounded  by   b.exp^-^c    }.      This    can  be    seen 
from  the    inequality 

|l^(2|cxsin(H)  |)exp[-|c2-x^)  | 

<   exp(-^     -X     +p|c|x) 


<   e 


12      3    2- 


(87) 


(2)   The  contribution  to  (86)  from  the  part  x  <  1  is 
1„2- 


bounded  by  b*exp[-2i<;  ).   This  follows  from 
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12       2 


12         2,1 


|l^(2|cxsln(H)  |)exp(-|c''-x^)  I    <    exp  (-^o^  -  x^  +i-|  c  | )     .       (88) 

Because    of   this,    we    can   divide   the   domain   of   Integration 
in    (86)    as  cS  and  .£'    where 

j&  =    {  |sin(H)  I    >   ^     and      x    >   l} 

and„^'  is  the  complement  of  ,3j  and  we  can  write  for  large  |c| 


r2(c,k)|   <   I  d3gdx{|gr^[|g|^+x^](|c|x)-^e-^} 


1.2- 


+  0[exp(-^^)] 


(89) 


Q  =  |c^  +  i-n^  +  x^  -  2x  I  csln  (g)  |  , 


(90) 


12^ 


where  the  ,JS  ' -contribution  is  written  as  0[exp(-TrC  )  ]  on 
account  of  (87)  and  (88),  and  we  have  used  large-c  asymptotic 


estimate  for  1   (2|cxsin  @|)  . 

It  is  a  simple  matter  to  obtain  from  (90)  the 
inequality 


Q   >   |c^cos^  +  ^^  +  (x-|csin(l{ 


which  when  substituted  into  (89)  yields 


(91) 


|r„(c,k)|  <  b(l+c2)*Y  1 


(92) 


Again,  we  can  use  (75)  and  partially  integrate  g,, 
and  the  result  is 


r„(c,k)|  <  b(l+c2)2Y-i  |1|  . 


(93) 


51 


Now  combining  (92)  and  (93)  yields  the  final  result 
for   rp(Cjk)  as 

|r„(c,k)|   <   bd+c^)"''  (1+k^)"''  .  (94) 


It  should  be  noted  that  for  differentiation  with 
respect  to  g.   we  have  to  consider  the  case  |sin  (h)  |>  -jt 
and  I  sin  (h)  |  <  -r-   separately.   For  |  sin  (e)  \    >  -r-  ,  we  note 
that  |cxg|   can  be  differentiated  with  respect  to  g  without 
difficulty;  the  differentiation  only  adds  an  extra  factor 
of  |c|  which  we  have  in  (93)'   On  the  other  hand,  for 
|sin  (h)  I  <  jj-  ,  we  use 

2    P   1/2 
I  8  lQ(2|cxg|cot  0)|<  b(l+c  cot  9)    lQ(2|cxg|cot  9).   (95) 

2    2   1/2 
Where  the  factor  (1+c  cot  9)     in  (95)  is  covered  by 

inequality  (87). 

Now  to  obtain  estimates  for  the  derivatives  of 

rp(c,k)  with  respect  to  c  and  k  we  again  have  to  consider 

I  sin  (h)  I  >_  -jT   and  |sin  (^  \    <_  ■^■   separately.  For  |sin  (h)  |^  |j-  , 

we  can  differentiate  |cxg|  directly  and  for  |sin  (h)  |  <   -^   . 

We  use 

|9^8glQ(2|cxg|cot  9)1  <  bd+c^+g^cot^g)^^!"^  Iq  ,      (96) 

where  the  factors  in  (96)  are  again  covered  by  (87). 
And  the   result  can  be  obtained  in  the  same  way  as  the 
procedure  for  the  estimate  (9^).   In  fact,  we  have 

|8^9JJr2(c,k)i   <   bd+c^)"''  (1+k^)"''  .  (97) 
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for  |m|  ^0,  |n|  >  0 


Summing  the  results  for  r, (c,k),  rp(c,k)  and  r^Ccjk) 


we  have 


|8j3j;r(c,k)|   <   b(l+c2)-^i+k2)-^  , 


(98) 


for  |m|  >0,  |n|  >0 


(C)   Estimates  for  C(c)q(c,k) 


Since  cCc)  =  0  for  |c|  >  2c  ,  we  only  have  to 
consider  q(Cjk)  for  |c|  <  2c  .   We  first  obtain  from  (17) 

Req(c,k)  >_  exp  (-ilc^-1)      d^g|g|'^    deg  (  G  ,a)slne  (secO  )  ^"^^ 


g 


<1 


2^.  2, 


r -,  ,   T   -g  tan  9 

[1  -  cosk'gj'e  ^ 


(99) 


We  now  substitute  (9)  into  (99)  and  note  that  the 
contribution  due  to  B(tan0)  In  (9)  gives  a  bounded  function 
b(c jk) J  I.e. 


Req(c,k)    >_  exp  (-i^c      -  1) 


d^g 


00 

r 


Islli 


3-2v         2 


X    [1-cosk.g]    +   b(c,k) 
k 


>   k^^exp (-4c^    -    1) 


-,      -l-2vrslny         ,"1 
dyy  |_-^  -  ij 


+   b(c,k) 


(100) 


or,    simply 


Req(c,k)       >      bk 


2v      1 


N 


(101) 
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for  |c|  <  2c   if  N  is  taken  large  enough. 

The  derivatives  of  q(c,k)  can  be  estimated  likewise 
and  the  result  is:   for  |m|  >_  0  we  have 


8"^q(c,k)|   <   b(l+k^)^ 


(102) 


and  for  |m|  >_  0   and   |n|  ^1,  we  have 


<     b(l+k^)~^ 


(103) 


where  both  (102)  and  (103)  are  valid  only  for  |c|  <  2c  . 


(D)   Estimates  for  v(Cjk) 

We  finally  come  to  the  most  important  term  v(c,k) 
The  procedure  we  use  is  somewhat  similar  to  that  for  the 
estimates  of  rp(c,k)  in  the  subsection  (B). 

We  first  write  (24)  out  explicitly  as 


v(c,k) 


,3  ^ n   ik'gs 
d-^gd-e   ^) 


dec(e)G(tane) 


X  [|gtane|^-|c|^]exp(-|c^-i-n^)  , 


(104) 


CO)  =  (sine)   (sece)^B(e,a)  , 


(105) 


where  G(tane)  is  given  in  (7). 

We  then  introduce  the  variable  x  =  gtanO  and  write 


C(6)d6  =  c(x,|g|)dx 


(106) 
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and  from  (9)  we  have  the  following  bound 

|iU.|g|)l      1     blgr^[l+    (x/|g|)2+2vj    .  (107) 

Substituting    (107)    into    (104)    yields 

00 

|v(c,k)|    <   b    I   d3g|l-e^^*S|    I  dx|grl[l   +    (x/|g|)2+2vj 

0 

[x^-|c|^)    I^(2|c|xsin(H))exp(-|c2-|fl2-x2)       (108) 

Again  we  observe  that  contribution  to  (108)  due 
to  the  parts  of  |sin@|  <  jr   and  x  <  1  are  bounded  by 
exp  (-ttC  J;  the  argument  is  the  same  as  those  given  in  (87) 
and  (88).   Consequently,  we  only  have  to  consider  the  domain 
(^   used  in  (B) : 

X'=  {|sin(D  I  >  ^  ,   X  >  1}  , 

where  for  large  c  the  function  I   can  be  asymptotically 
expanded;  the  result  is 

|v(c,k)|  <  b  I  d3gdx/|l-e^^'S||g|-iri+  (x/|g|)^'^2^"| 


X  |(x^-|cr)|(|c|x)-^e-^} 
+  o[k2^exp(-lc2)] 


(109) 


where  Q  satisfies 


Q  =  |c^.  +  |-n^  +  x^  -  2x  I  csin  (Ti 


1  2 


1  2 


>  jc^cos'^Ch)  +  ^'^  +  (x-  |csin(H)  I  j'^       (110) 
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Now  for  large  |c|,  we  observe  that  Q  Is  large  unless 
=  0,   sin  (jy  =  1   and  x  =  |c|.   We  therefore  collect 
a  factor  |c|~   from  the  integration  sin  (H) d  (H) ,  a  |c|~ 
factor  from  (jclx)""^,  a  factor  |  c  |  ^~   from  |  x^  -  |c|^|  and 
a  factor  Ic]^"^^^  from  x^"*"^^.   The  net  result  for  (109)  is 

|v(c,k)i   <   b(l+c^)~''[k2+l]^  (111) 

which  is  the  desired  result.   We  can  make  the  above  computa- 
tion more  rigorous  by  using  y  =  |c|cos(2)j   t  =  x  «•  |  c  sin  (^  | 
as  new  integration  variables  in  (109)  and  expand  for  large  c; 
the  details  are  straightforward  but  tedious  and  we  shall  not 
reproduce  them  here . 

To  obtain  estimates  for  the  derivatives  of  v(c,k)  with 
respect  to  c,  we  again  have  to  consider  the  cases  |  sin  (^  |  <  jj- 

1  "I 

and  I  sin  (H)  I  >  jj-  separately .   For  |  sin  (H)  |  <  jj-,  we  use 


9™I  (2  I  cxsin  (S)  I )  I   <   b(l+c^+x^)'^I  (2  |  cxsin  (h)  |  )     (112) 


and  the  factor  (1  +  c  +x  )   is  covered  by  the  factor  exp  (-^  ) 
in  (109). 

For  I  sin  (H)  I  >  jj-,  the  procedure  goes  as  follows. 
First  consider  the  term 

Z  =  I  (2  I  cxsin  (H)  |)exp(-p-c  -^^  -x  ) 


2tt 
2tt 


e  "d(})  (113) 
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in  eqn.  (108)  where 


n  =  hi^    +   |ri^  +  x^  -  2cx  sln(H)cos(|) 


=  (c  cos  @  +  |g)^  +  ^^  +  (x-|csln(H)  |cos({))^ 

+  c^sin^(H)sln2<}.  ^  (ll4) 

We  then  differentiate  with  respect  to  c  and  obtain 

||3c"l    i   I  (ccos®   +|g)|   +  I  (x-|csin(H)  |  cos(())  1 

+    |csin(H)  sin^l  (115) 

which   shows    that 

2tt 


I'c^l      <-     ^1     ^-^"^» 


1_2    1_2    1,2- 


<   bI^(|cxsin(H)  |)exp(-^'^-^n    -^x^)     .         (Il6) 


On  the  other  hand,  since  |sin(H)  |  >  j^-,    repeated 
differentiations,  which  give  factors  like  |sin@|~  ,  cause 
no  problem.   And  a  similar  argument  as  that  given  above 
shows  that 


|3'"z|   1  bl  (Icxsin®  |)exp(-^2-^n^-|-''^)  •        (117) 


Now  differentiate  v(c,k)  in  (104)  and  using  (117), 
we  obtain  the  same  type  of  integrals  as  (108);  in  fact  we 
would  have,  exactly  the  same  integral  if  we  rescale  c  into 
/2  c,  n  into  /2  n  and  x  into  /2  x  .   We  therefore  arrive  at 
the  same  estimate  as  given  in  (111),  namely 
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8%(c,k)|   <   b(l+k2)'^^(l+c^)  ^  (118) 


for  |m|  >  0. 

Finally,  we  consider  differentiation  of  v(c,k)  with 
respect  to  k.   Since  differentiation  with  respect  to  k.  only 
brings  down  a  factor  g.  from  (1-e   ^),  the  same  procedure 
as  In  deriving  (ll8)  goes  through  without  change  and  we  have 

|3^3;jv(c,k)|   <   bCl+c^)"''  (119) 

for  |n|  >^  1  where  we  do  not  have  the  factor  k  ^  because  of  the 
presence  of  the  extra  factor  g.. 

But  on  the  other  hand,  since  the  k-dif ferentiatlon 
brings  an  extra  factor  g.,  we  can  then  use  (75)  and  integrate 
by  parts  with  respect  to  g,,  .   In  doing  this,  we  shall  raise 
the  power  of  |c|  by  1,  which  comes  from  differentiating  ^ 
in  (114)  with  g,,  ,  and  the  result  is 

|9^3;jv(c,k)|   <   b|krl(l  +  c2)-^+^(l+k2)^       (120) 

where  |n|  ^1.   Now  we  combine  (119)  and  (120)  to  obtain  the 
desired  result 

|3";3%(c,k)|   <   bd+k^)   ^d  +  c^)   ^  (121) 


c  k 


(1-2V)  /TOON 


for   n   >  1, 
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(E)  Estimates  for  s„(c,k) 


We  write  (3^)  as 


s„ (c,k)  =  |c 


*Tr  3 

B(e,a)(sin9)^^(2t)^exp[-|-tcos^ec^) 


X  h  -  exp(-|-tk^cos^9) Ide  +  b(c,k)  (123) 

where  the  function  b(c,k)  and  its  derivatives  with  respect  to 

c  and  k  are  bounded  functions  for  all  c  and  k.   If  we  use 

2    2 
tk  cos  6  as  a  new  integration  variable  for  the  Integral  in 

(123),  we  find  that 

Resj^(c,k)  >  b|c|^k^^exp(  -  c^/k^}  -  ^N  (124) 

if  we  take  N  large  enough. 

On  the  other  hand,  differentiations  with  c  and  k 

2 
bring  down  extra  factor  of  cos  9,  and  consequently  we  have 


9^s^(c,k)|  <  b  |m|  >  1  (125) 


9^8;js^(c,k)|  <  bd+k^)"''    |n|  >1  (126) 


(F)  Estimates  for  c(c)s(c,k) 
First  we  can  give  a  bound  for  s(c,k)  by  writing  (30) 


as 


lsCc,k)|  1  2|c|^r  B»(G)A  dOe 

(127) 
A  =  [1  -  J  (t  sin  29|cxk| )] 

+|{1  -  exp[-  |t  Gos^e(k^+2ic -k)]} 1 
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We  can  use  |cxk|cos  9  as  the  new  integration  variable 
In  (127)  for  the  first  term  of  A  and  cos  6|k|  for  the  second 
bracketed  term  In  A.   Since  the  singularity  of  6*  (6)  at  6  =  ^tt 
Is  given  by  (cos  0)~  ~  ^  ,  we  obtain  from  (126): 

?(c)|s(c,k)|  <  b|k|^^  (128) 

This  shows  that  If  e  In  (37)  Is  small  enough,  then 

C(c)  Re[q(c,k)  -  es(c,k)]  ^|b|k|^^  "  S"  ^  (129) 

which  follows  from  (101). 

Now  for  the  derivatives  of  s(c,k)  with  respect  to  c  and  k, 
we  use  the  relation 

1(3^)"^  J^(x)l  =  |x-^  J^(x)|  ib  (130) 

where  x  denotes  t. sin  29|cxk| .   We  note  that  the  derivatives  of 

2 
X   with  respect  to  c  and  k  can  be  easily  calculated.   The  end 

result  may  be  stated  as 

|3^[C(c)s(c,k)]|  <  bd+k^)"",  (131) 

|3^  9jJ[c(c)s(c,k)]|  ^bd+k^)""".  (132) 

(G)   Summary 

We  now  use  (37)  and  rewrite  the  symbol  a£(c,k)  as 

a^Cc,k)  =  es^(c,k)  +  c(c)[q(c,k)  -  es(c,k)] 

+  [1  -  ^(c)][(l-e)s^(c,k)  +  v(c,k)] 
+  r(c,k)  +  u(c)  (133) 

The  results  in  the  preceding  subsections  (A)  -  (F)  may  be 

summarized  as  follows. 

2     2 

For  c   <  k  ,  we  obtain  from  (12^) 

Re[(l-e)s^(c,k)]  >  b  |  c  |  ^  |  k  |  ^"^  -  g-N  (13^) 
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2  2 

and    for   c      >    k    ,    inequality    C65)    Indicates   that 

uCc)    >    |blcnk|2v    ^   ib|cr+2v    _    1   ^  (^35) 

which  together   with    (13^)    yield 

Re[Cl-e)Sjj^(c,k)    +   |uCc)]    >    b(c^+k^)^   -   g^  N  (136) 

for  all  c  and  k.   We  next  consider  (111)  and  note  that  if  lc| 
is  large  enough,  then 

Re[(l-e)s^(c,k)  +  |u(c)  +  v(c,k)]  >  bCc^+k^)"^  -  Jn  (137) 

and  we  shall  take  c   to  be  so  large  that  (137)  is  valid  for 

IcI  >  c  .   We  note  that  the  choice  of  c   can  be  made  independent 
'  '     o  .0 


of  I 


On  the  other  hand,  for  |c|  <  2c  ,  we  can  use  (129).   The  end 


result  is 

Re  aj^(c,k)  >  Re[es  (c,k)  +  |-u(c)] 

2Y 

+  N  +  b(l+c^)   (c^+k^)^  (138) 

for  all  c  and  k. 

We  then  note  that  by  using  the  formula  of  Hilb  for  P„(sin  G) 
in  s.(c,k),  we  can  show  that 

Re[es^(c,k)  +  |u(c)]  >  h{l  +  '^)^^    -  ^N  (139) 

for  all  c  and  k,  and  its  proof  is  similar  to  the  arguments  used 
in  Sec.  5  of  Part  I.   Combination  of  (138)  and  (139)  yields  (^2). 

Finally  it  is  noted  that  (43)  follows  from  (70),  (98),  (102), 
(118),  (125)  and  (131)  while  (44)  follows  from  (98),  (103),  (121), 
(126)  and  (132)  . 
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5.   The  Eigenvalue  Problem 

With  the  results  in  (42),  (43)  and  (44)  established,  we 
are  now  In  a  position  to  Inverst  the  operator  Aj  in  (40).   In 
fact  we  can  introduce 


e^^*^[a^(c,k)]-^$(k)d3k 


(140) 


and  show  that 


BA„  =  1  +  R  , 


(141) 


IIRII   <   b/N  , 


(142) 


where  II  Rll  is  the  o^  -norm  of  R,  and  b  independent  of  £. 

The  proof  of  (l4l)  and  (l42)  is  exactly  the  same  as 
that  given  in  Sec.  6  of  Part  1.   From  this  result,  it  follows 
then 

A-1  =  (BA^)-^B  (143) 

—  1  -2      2 

which  implies  that  A^   is  a  compact  operator  f rom  ^.  to -5^  , 

-2     -2 
since  B  is  a  compact  operator  f  rom  ,?£  to  'jt     according  to 

Theorem  2  in  Sec.  3  of  Part  I. 

From  (42)  we  also  note  that 

I [a^(c,k)]~^|  ^  0 

as  £  ->  00  and  consequently  we  have,  as  £  ->  0°, 

-111 


A 


I 


0 


(144) 


We  therefore  arrive  at  the  following  conclusions. 
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(1)  A   is  a  compact  operator  from  jiL  to  ^  . 

(2)  L~   is  a  compact  operator  from  ^  to  ^  ,  where  ^ 
Is  the  subspace  of  X     which  is  orthogonal  to  the  five  eigen- 
functions  of  L  for  the  eigenvalue  zero  (i.e.  orthogonal  to 
the  five  collislonal  invariants  with  the  proper  weight  factor) 

(3)  The  spectrum  of  L  is  discrete. 

2 

(4)  The  eigenf unctions  of  L  are  all  )t  -functions  and 

2 
they  are  complete  in  ^  since  L  is  symmetric  and  non-negative 

and  therefore  essentially  self-adjoint. 

(5)  The  X  -eigenf unctions  of  L  in  fact  belong  to  ^; 
this  will  be  proved  in  the  next  section. 
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6.   Properties  of  the  Eigenf unctions 

We  shall  now  show  that  the  eigenfunctions  <}>  of  L  in 
fact  belong  toj/;  I.e.  they  are  C  -functions  having  a  decay 
at  Infinity  faster  than  any  inverse  powers  of  c. 

We  first  show  that  the  eigenfunctions  (J)  are  C  -functions 
To  this  end,  we  recall  eqn.  (40)  and  write 

A(J)   =   Aj^cf)   =  Y^    ,  (145) 


Yi|^   = 


;^^*'=a,(c,k)(l+k2)-^^$(k)dk  ,       (146) 


$(k)   =   (l+k2)^^$(k)   .  (147) 

Then  we  can  invert  Y  in  the  same  way  as  we  invert  A. 

—  1  2      2 

and  Y   is  a  bounded  operator  from^  to  56  .   Consequently, 

since  ()>  e  ^^  , 

shows    that    ij;   e  ^^    or    (l+k^)^^$(k)    e^^. 
Next   we   write 

rp      =      X<p      =       I    e^^*^(l+k2)^^$(k)dk  (148) 

and   apply   the   operator   X   to    (l45).      This    gives 

X\p      =      YXip   +    [X,Y]i|;  (149) 

where  [X,Y]  is  the  commutator  of  X  and  Y. 

Now  applying  Y~  to  the  last  equation,  we  obtain 

Xip   =   AY~^ip  -  Y~^lX,Y^\p  (150) 
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It  can  be  shown  that  the  operator  Y   [X,y]  Is  a 
bounded  operator  from^     to^  ;  the  proof  is  similar  to  the 
procedure  in  Sec.  6  of  Part  I  using  the  reverse  operator 
and  multiple  symbols  and  we  shall  not  give  it  here.   The 
conclusion  is  then  X-^i   e^     or  (l+k2)^$(k)  e  5^.  . 

Proceeding  in  this  manner  with  mathematical  Induction 
leads  to  the  result  that 

(l+k2)^'^^$(k)  e^2  ^^^^^ 

for  all  positive  integer  n.   This  shows  that  ())(c)  is  a 

CO 

C  -function. 

In  this  proof,  the  main  property  we  use  is  the 
growth  (1+k  )   in  k  in  a.Ccjk)  in  (42).   Since  a„(c,k)  also 
has  a  growth  (1+c  )   in  c ,  we  can  use  similar  arguments 
and  show  that 

(l+c2)^'^^(f(c)  £56^  ,  (152) 

but  the  proof  is  slightly  more  complicated  than  that  for  (151 ) 
and  will  not  be  reproduced  here. 

It  is  further  mentioned  if  we  use  a  mixture  of  the 
arguments  used  in  establishing  (151)  and  (152)  we  can  show  that 
(1+c  )   ^^    h  ^   belongs  to  X  ;  i.e.  (f)  e  ^  » 
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7.   Concluding  Remarks 

Finally  we  make  a  remark  concerning  the  part  of  the 
collision  Integral  conventionally  written  as 


C^[f]  =   Id^c^  exp(-c^) 


de  dQ  a(|g| ,e)[f^-f^] 
1  -2. 


(153) 


If  we  use  ({i(c)  =  f(c)  exp(-p-  c  )  then  the  operator  L,  , 

12  12 

L-^(i,    =    exp(-2  c  )C-|^[(j)(c)  expC^c  )]  , 

Is  a  compact  operator  from  X     to  ^  ,  which  can  be  seen  from 
the  estimates  in  (98)  and  Theorem  3  In  Sec.  3  of  Part  I. 
This  shows  that  all  the  difficulty  in  the  Boltzmann  Collision 
operator  stems  from  the  term  involving  [f'-f]. 
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